Entanglement Entropy in Loop Quantum Gravity 
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The entanglement entropy between quantum fields inside and outside a black hole horizon is a 
promising candidate for the microscopic origin of black hole entropy. We show that the entanglement 
entropy may be defined in loop quantum gravity, and compute its value for spin network states. 
The entanglement entropy for an arbitrary region of space is expressed as a sum over punctures 
where the spin network intersects the region's boundary. Our result agrees asymptotically with 
results previously obtained from the isolated horizon framework, and we give a justification for this 
agreement. We conclude by proposing a new method for studying corrections to the area law and 
its implications for quantum corrections to the gravitational action. 
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INTRODUCTION 

Quantum field theory in curved space indicates that 
black holes radiate thermally with an entropy given by 
the Bekenstein-Hawking formula 
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where A is the area of the black hole event horizon [l|, 
2] . A similar entropy is present for acceleration horizons 
in flat space J3(] and cosmological event horizons in dc 
Sitter space [4|]. It has therefore been suggested that 
every causal horizon possesses an entropy proportional 
to its area, so that ([I]) represents a universal horizon 
entropy density [Bj]. 

Although there is no universally agreed-upon source of 
black hole entropy, a strong candidate is the entangle- 
ment entropy of quantum fields inside and outside the 
black hole horizon Q. This approach is based on the 
observation that entanglement between fields inside and 
outside the horizon causes any globally pure state to be- 
come mixed when restricted to the exterior of the black 
hole. The entanglement entropy has the main features 
expected from a statistical description of black hole en- 
tropy: it scales like the area of the horizon in the presence 
of a Planck-scale cutoff 0, 0] and under the assumption 
of causality it satisfies the generalized second law of ther- 
modynamics 0. 

In the previously studied case of a scalar field on a flat 
background, the entanglement entropy of a spherical re- 
gion diverges in the absence of an ultraviolet cutoff. It 
has been suggested that quantum gravity could act as an 
ultraviolet cutoff, rendering the entanglement entropy fi- 
nite. Moreover, if the gravitational field is quantized it 
will also contribute to the entanglement entropy. There- 
fore any investigation of the relation between entangle- 
ment entropy and the Bekenstein-Hawking entropy must 
be done in a theory of quantum gravity. 

The goal of this paper is to study the entanglement 



entropy of the gravitational field within the framework of 
loop quantum gravity. The state of the gravitational field 
will be described by a spin network state 9, 10]. The spin 
network states are chosen because of their interpretation 
as states of discrete geometry. In particular, spin network 
states are eigenstates of the area operator, allowing a 
comparison between the entanglement entropy and the 
horizon area. 

Entanglement entropy has previously been computed 
in loop quantum gravity for black hole coherent states 
in spherically symmetric spacetimes with apparent hori- 
In contrast with previous work the present 
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result is obtained without assumption of symmetry or of 
particular boundary conditions at the horizon. 



ENTANGLEMENT ENTROPY AND THE 
SCHMIDT DECOMPOSITION 

We first define the entanglement entropy and review 
some of its well-known properties, including its relation 
to the Schmidt decomposition. 

Let M be a spacetime manifold decomposed as M. = 
R x £ into time and space. Let denote the space 
of wave functionals ip : <!>£ — > C where $s is a suitably 
defined space of field configurations. A partition of space 
£ = ft U Q gives a tensor product decomposition of the 
Hilbcrt space 



(2) 



Let \tp) € Ht, be a pure state. For each region ft there 
is an associated mixed state in 7ir> given by tracing over 
degrees of freedom in f2, 



(3) 



This state encodes all information that can be obtained 
about \ip) by performing measurements localised in f2. 
The entanglement entropy of the region Q is defined as 



S E (p) = S(p(Sl)) 



(4) 
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where S(p) = — Tr(plogp) is the von Neumann entropy. 

Definition (The Schmidt Decomposition). Let £ 
Tin <8> Hq- Then there exist orthonormal sets {\rp^ )} C 

Tin and ||v4 2 ^} c an d positive real numbers {A,} 
such that 



(5) 



The numbers |\/A7} are called the Schmidt coefficients, 
and the number of terms in the sum is the Schmidt rank. 

From the Schmidt decomposition, we can compute the 
diagonal form of the reduced density matrices 
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This shows that the two reduced density matrices have 
the same nonzero spectrum. It follows, as is well known, 
that the entanglement entropy is symmetric and can be 
computed from the Schmidt coefficients 



(8) 
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In order to compute the entanglement entropy of a spin 
network state, it is therefore sufficient to compute its 
Schmidt decomposition. 



• rh — (mi, . . . , tun) where m n is a vector in the spin 
J n representation space Vj n = C 2J " +1 . 

Let A £ A be a generalized connection, A (7) the holon- 
omy along curve 7 and R J (^(7)) its spin j representation. 
Then the extended spin network state IS*) is defined by 

(A\S) = ((g)^(A(7^))®(g)m n ) o (® in ) (9) 

\e=l 71=1 / \n=l / 

Throughout we will assume that all intertwiners i n are 
normalized so that Tr(i* i n ) = 1 and vectors m n are nor- 
malized so that ||m n || = \/2J n + 1. 

An important property of these states that we will use 
is their orthogonality. Suppose that S — (T,j,J,i, rh) 
and S' = (r, j, J, 1, rh') so that S and S' differ only by 
their set of vectors. Then 

X (m n \m' n ) 



(s\s') = n 
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The second important property of the extended spin 
network states we will use is the freedom to insert ver- 
tices. Let S be an extended spin network. Given an arbi- 
trary point v on the curve jl, we can split 7^ into = 
l'h 7l so tnat 7l an d 7l mee t at v - Let T' be the graph 
with links 71, ... , jl-i, IliIl an( i vertices v\, . . . , vn, v. 
Let S' be the extended spin network (T' ,j' , J' ,1' ,fh') 
where / = (ji, . . . , j L , j L ), J' = ( Ji, . . . , J„, 0), 1' = 



mi, 



, m„, 1). Then the re- 



sulting state 1 5") is equivalent to \S), \S'} = \S). 



Spin network states 



ENTANGLEMENT ENTROPY OF SPIN 
NETWORK STATES 

In loop quantum gravity, the space of fields is the space 
A of generalized connections on E. We will consider 
the Hilbert space of cylindrical functions, which is 
spanned by the extended spin network states [9(. 

An extended spin network is a tuple S = (T, 3, J, i, rh) 
where 



• r is a graph in E consisting of nodes V\ , . 
links 71, . . . ,7 L , 



, vn and 



• J = (h, ■ ■ -Jl) where j e £ {\, 1, §, . . .} labels a 
non-trivial irreducible representation of SU(2), 

• J = (Ji, . . . , Jn) where Jn £ {0, -, 1, |, . . .} la- 
bels a possibly trivial irreducible representation of 
SU(2), 

• 1= . . . , ijv) where i n is an intertwining operator 
from the representations of all incoming edges and 
the spin J n representation to the representations of 
all outgoing edges, and 



We now consider the entanglement entropy of a general 
spin network state. A spin network is an extended spin 
network for which the representation attached to each 
node is trivial, in other words J n = 0, m n = 1. 

Suppose that f2 is a subset of E whose boundary dfl 
intersects T only at links. We can insert nodes with in- 
tertwiners at all points where the graph T intersects the 
boundary dQ. Let P be the number of points where T in- 
tersects <9f2, Nq the number of vertices of T in fi, and Nq 
the number of vertices in fi. We can partition the nodes 
vi, . . . ,vn so that v n £ dfl for n = 1, . . . , P, v n € for 
n = P+l, . . . , P+Nq, v n £ n for n = P+N Q + 1, ...,N. 
Similarly, the links 71, . . . ,7z, can be partitioned so that 

€ fl for £ = 1, . . . , Ln and 7^ € f2 for £ = Lq+1, . . . , L. 

By the construction for inserting vertices, S has i p — 
—j==I for p — 1, . . . , P where j p denotes the spin of 

the edges incident to node p. Each of these intertwin- 



ers i p may be expanded in an orthogonal basis of Vj, 
{ei,...,e 2 j p+ i}. 



2Jp + l 

■ E 

a„ = l 



(11) 
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We can apply this decomposition to each of the inter- 
twiners on the boundary. Letting a — (ax, . . . , ap), 



(12) 



where the sum over a means to sum over all n-tuples 
(ax, . . . , a P ) with a p = 1, . . . , 2j p + 1. 

The state \Sq,o) is an extended spin network state 
with graph Tq consisting of the links j£, I = 1, . . . , Lq 
and vertices v n ,n — 1, . . . , P + Nq. The labels on the 
links are unchanged, j' e = ji for £ = 1,...,Lq. The 
existing nodes are also unchanged J' n — J n , i' n — i n and 
m' n = m n for n = P + 1, . . . , P + Nq . The inserted 
vertices on the boundary have J' — j p , i' = — . I 



and mi 



( , y/2j p + le ap for p = 1, . . . , P. 
By the orthogonality relation (jTOJ) , these states form 
an orthonormal set. Therefore (fT2|) is a Schmidt decom- 
position of |5). The Schmidt rank is 



Here the superscript a and the subscript b index the ma- 
trix elements of the intertwiner. Because the set of states 
Sq j and Sq, b^j are orthogonal, this formula reduces 
to the entanglement entropy of the tensor product of all 
boundary intertwiners. 

By additivity of the von Neumann entropy across ten- 
sor products, the entanglement entropy is 



P =i 



(17) 



In this equation i p is viewed as an entangled state be- 
tween the representations of edges incident from £1 and 
those incident from Q. 

Note that in the case where i p is a multiple of the 
identity this formula reduces to 5b (z p ) = log(2J p + 1) as 
in equation (fT4|) . 



RELATION TO THE ISOLATED HORIZON 
FRAMEWORK 



N=Y[(2j p + l) (13) 
p=i 

From (|8|) the entanglement entropy of |5) is 

p 

5 £ (r!)=^log(2j p + l) (14) 
P =i 

The Schmidt decomposition also allows us to compute 
the reduced density matrix corresponding to the region 

n. 



p(Q) 



— ^2\Sn,a){Sn,a\ 



(15) 



Note that although the individual states \Sq, a) trans- 
form non-trivially under a gauge transformation, the lin- 
ear combination (|15|) is gauge- invariant. 



Intertwiner on the boundary 

So far it has been assumed that no nodes of the spin 
network 5 lie on the boundary of f2, in which case the 
entanglement entropy depends only the edges of 5 that 
intersect the boundary. In the case where a node lies on 
the boundary, the entanglement entropy depends on the 
intertwiner assigned to the boundary node. 

In this case we can let P be the number of intertwiners 
on the boundary 



1 5) = 



E 

d.b 
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Sq, a^}® 


\p=i / g 
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The expression ([13]) for the Schmidt rank of the spin 
network state coincides with the dimension of the bound- 
ary Hilbert spaces defined in the isolated horizon frame- 
work 12 , 13| • The analogy between the isolated horizon 
framework and entanglement entropy was first suggested 
by Husain, who showed that the calculation of black hole 
entropy in the isolated horizon framework does not de- 
pend on the details of the boundary conditions 14j. To 



elucidate the relationship between these two theories it 
is necessary to first introduce the relevant aspects of the 
isolated horizon framework. 

In the isolated horizon framework, the horizon is 
treated as an inner boundary of space obeying the iso- 
lated horizon boundary conditions. The gravitational ac- 
tion acquires a surface term proportional to the action of 
Chern-Simons theory. The states are therefore elements 
of Tin <S> Ti-do, where Tin is the "bulk space" of cylindrical 
functions on fl and Hon is the "boundary space" of U(l) 
Chern-Simons theory on dQ. 

The isolated horizon boundary conditions are imple- 
mented as an operator equation on TLq ® Tlgsi restricting 
states to a subspace 



(18) 



Where V runs over all finite set of punctures labelled by 
spins j p . The space Hq is spanned by the open spin 
network states |5n,a), where Sq runs over all extended 
spin networks intersecting dfl in the set of labelled punc- 
tures V . H-qq is the space of Chern-Simons states on the 
punctured surface dil — V . 

A partial trace over the bulk space TLq is performed, 
yielding a maximally mixed state on Ti-g^- The entangle- 
ment entropy is therefore given by log dim "Hqq , where in 
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the limit of a large number of punctures logdim7ig f2 is 
given by 

p 

dimW^~n( 2 ^ + 1 ) ( 19 ) 

P =i 

which is the same as the Schmidt rank of a spin network 
that intersects dCl at the points V . 

To see why these quantities should agree, consider an 
arbitrary boundary theory with Hilbert space Tin^Tign- 
We will make the assumption that the boundary condi- 
tions should not affect physics in f2. Therefore for each 
spin network state \S) G Tin <8> Ti^ there should exist a 
state |S") G Tin €5 Tien that describes the same physics 
on the exterior 

TY W „(|SXS|)=TWIWI) ( 2 °) 

Now consider the mixed state of the boundary pgn — 
Tr-^ n (|S")(S"|). By equation ([3]), pan has the same 
nonzero spectrum as pn ■ Therefore the range of pan is a 
subspace whose dimension is given by the Schmidt rank 
(TT3"]) . Thus equation (fT^)) is a consequence of the fact that 
the rank of pn is the same as the rank of pgn- 

By imposing the requirement that the boundary con- 
ditions should not restrict the exterior Hilbert space, we 
have derived the relationship between the Schmidt rank 
of a spin network state and the dimension of the bound- 
ary Hilbert space. However equation p^|) only holds 
asymptotically; the dimension of Tig^ is less than the 
Schmidt rank. 

To understand the reason for disagreement, we note 
that a basis of the space Tig n is labelled by se- 
quences of half-integers x satisfying m p £ 
{— J p , —j p + 1, . . . , j p } and the additional spin projection 
constraint 

p 

^m p = (21) 

P =i 

Therefore the difference between the dimension and the 
Schmidt rank is entirely due to the spin projection con- 
straint. 

OUTLOOK 

We have computed the entanglement entropy of the 
gravitational field in loop quantum gravity for an arbi- 
trary region of space CI and spin network state \S). The 
entanglement entropy is a finite and extensive quantity 
that depends linearly on the number of punctures of the 
horizon. No assumptions have been made about the re- 
gion CI, so our result applies to all causal horizons. 

We have found an interesting relation between the en- 
tanglement entropy and the isolated horizon framework. 



What remains to be understood is the mismatch between 
the two theories due to the spin projection constraint. It 
would therefore be of interest to see whether this con- 
straint is necessary, or whether it can be relaxed by a 
less restrictive choice of boundary conditions. Of partic- 
ular interest is the relation between the spin projection 
constraint and the requirement that dCl have spherical 
topology. 

Finally, the fact that the entropy is extensive over the 
horizon allows thermodynamics to be applied locally to 
horizons in loop quantum gravity 1161 1. It has been 
shown that under the assumption that the horizon en- 
tropy is extensive, properties of the effective action for 
gravity can be computed from the entropy density. How- 
ever, in order to apply these arguments it is necessary to 
express the entropy as a function of geometric variables 
instead of spin networks. This suggests the question of 
whether there exists a geometric quantity Q such that 




Results in the context of classical general relativity sug- 
gest that Q is related to a Noether charge of the gravi- 
tational effective action Together, these results sug- 
gest that horizon thermodynamics could be a powerful 
tool for studying quantum corrections to the gravita- 
tional action in loop quantum gravity. 
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